Coupled Dynamics of Voltage and Calcium in Paced Cardiac Cells 
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We investigate numerically and analytically the coupled dynamics of transmembrane voltage 
and intracellular calcium cycling in paced cardiac cells using a detailed physiological model and 
its reduction to a three-dimensional discrete map. The results provide a theoretical framework 
to interpret various experimentally observed modes of instability ranging from electromechanically 
concordant and discordant alternans to quasiperiodic oscillations of voltage and calcium. 



Over the last decade, there has been a growing recog- 
nition that dynamic instability of the cardiac action po- 
tential can play a crucial role in the initiation of life- 
threatening arrhythmias Q, 0, IS 0, IS @ • Most studies to 
date have focused on the dynamics of the transmembrane 
voltage governed by the standard equation 
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where C m is the membrane capacitance, 7i on is the total 
membrane current, which is the sum of the individual 
currents for Na + , K + , and Ca 2+ ions depicted schemati- 
cally in Fig. ^ an d -^cxt is the external current represent- 
ing a sequence of suprathreshold stimuli equally spaced 
in time by T, the pacing period. A widely used approach 
to model the nonlinear dynamics of voltage is the one- 
dimensional discrete map A n+ \ = f(T — A n ) which re- 
lates the action potential duration (APD) at two sub- 
sequent beats via the restitution curve, A n+ i = f(D n ), 
where D n is the interval between the end of the previ- 
ous action potential and the next 0, 0, IS EL IS El ■ The 
periodic fixed point of this map corresponding to the sta- 
ble 1:1 rhythm undergoes a period-doubling instability to 
alternans, a sequence LSLS... of long (L) and short (S) 
APD, when the slope of the restitution curve is > 1. 

Even though this map has been successful to model the 
unstable dynamics of voltage in some ionic models 0] and 
experiments its predictions are inconsistent with a 
wide range of observations [a, Ifl □> IS| • F° r example, Hall 
et al. ; f| found that alternans can be absent even when 
the slope of the restitution curve is significantly larger 
than one, and conversely alternans are observed under 
ischemic conditions in which the restitution curve is flat 
. An important limitation of the restitution relation- 
ship is highlighted by recent experimental 0, IS El an d 
theoretical studies |ll| which suggest that alternans may 
result from an instability of intracellular calcium cycling. 
The coupled nonlinear dynamics of voltage and calcium, 
however, remains largely unexplored. 

In this letter, we investigate this dynamics by a nu- 
merical study of a detailed physiological model and an 
analysis of the dynamics based on iterated maps. The 
model consists of Eq. ^ with membrane currents (Fig. 
^) modeled based on modifications by Fox et al. 12] of 
the Luo-Rudy currents |l3j] , coupled to equations from a 
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FIG. 1: Illustration of currents that control the dynamics of 
voltage and intracellular calcium cycling. 



recent model of calcium cycling |ll| 
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where c, 



and ca arc the concentrations of free Ca 2+ in 



a thin layer just below the cell membrane (submem- 
brane space), in the bulk myoplasm, and the sarcoplas- 
mic recticulum (SR), with volumes v 8 , V{, and v aT , re- 
spectively, where the SR volume includes both the junc- 
tional SR (JSR) and the network SR (NSR); d } is the 
average JSR concentration in the whole cell as defined 
in Ref. |ll|. The concentrations c s and c$ are in units 
of fiM, whereas Cj and c'- are in units of [iMv ST /Vi. All 
Ca 2+ fluxes are divided by Vi and have units of /iM/s. 
Instantaneous buffering of calcium to SR and calmod- 
ulin sites in Vi and v s is accounted for by the functions 
P s = 0(c s ) and /3i = P(c{), and the currents 1^' describe 
time-dependent buffering to troponin C |Tlj . 
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Calcium release from the SR is triggered by calcium 
entry into the cell via calcium-induced-calcium-release 
(CICR) Release occurs at a very large number of 

junctions where several L-type Ca channels (lea) an d 
a few release channels (ryanodine receptors; RyRs) face 
each other in close proximity. Only one of these junctions 
is shown in Fig. ^for clarity. The total release current for 

the whole cell is the sum I Ie \ — Y^k=i ^rel' 01 l° ca l cur- 
rents I* ei at each junction where release channels are acti- 
vated. Active junctions appear as bright localized spots, 
or "sparks" , in confocal microscope imaging of calcium 
activity |ljj. The number of sparks N(t) varies in time 
since sparks are recruited stochastically and extinguish. 
The model takes into account this spatially localized na- 
ture of release and the dynamical equation for the release 
current (Eq. |SJ) captures phenomenologically three key 
experimental observations: (i) sparks are recruited at a 
rate proportional to the whole cell lea, or N ~ lea [Tfij |- 
which insures that calcium release is graded with respect 
to calcium entry 0, 0, (ii) the spark life-time r r is 
approximately constant, and (iii) the amount of calcium 
released increases with SR concentration (SR-load) [Ts|. 

Instability mechanisms. Ca 2+ alternans, a period- 
doubling sequence Isls... of large (I) and small (s) cal- 
cium transient (Isls peak c,-), can occur independently of 
voltage alternans in experiments with a single cell paced 
with a periodic voltage waveform 0. Both theoretical 
analyses |0, and recent experiments support that 
a steep dependence of release on SR-load is the under- 
lying mechanism of these alternans. The sensitivity of 
release to SR-load is controlled in the model by the slope 
of the function Q(c'j) at high load 

u = dQ/dc' 3 . (7) 

For a large enough slope, the model produces Ca 2+ alter- 
nans when paced with a periodic voltage waveform [Tl| 
as in the experiments of Ref. . 

Steep APD-restitution in the absence of Ca 2+ alter- 
nans can also induce APD alternans. This steepness is 
especially sensitive to the recovery from inactivation of 
the calcium current \VA Il3| 

^Ca = d f / Ca «Ca, (8) 

where ica is the single channel current and d (/) is a fast 
(slow) voltage-dependent activation (inactivation) gate. 
For the intermediate range of pacing rates studied in the 
present work, increasing the time constant rt of the / 
gate in the equation / = (/oo(^) — f )/ T f steepens APD- 
restitution and promotes voltage alternans. 

Voltage- calcium coupling. The mutual influence of 
voltage and calcium during the action potential is con- 
trolled by the membrane currents that depend on intra- 
cellular calcium concentration. These include ica an d 
the sodium-calcium exchanger iNaCa- A crucial property 
is that a change in the magnitude of the calcium tran- 
sient has opposite effects on these currents with respect 
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FIG. 2: Illustration of the effect of an increase in the magni- 
tude of the calcium transient, which can prolong or shorten 
the APD for (a) positive and (b) negative coupling, respec- 
tively. The sign of the coupling depends on the relative con- 
tributions of Jca and ^NaCa to the APD. The solid or dashed 
lines correspond to the same beat. 



to prolonging or shortening the APD. A larger calcium 
transient following a larger release enhances inactivation 
of lea v i a the calcium-dependent gate /c a , a nd hence 
shortens the APD, but increases the chemical driving 
force for Ca 2+ extrusion from the cell via the exchanger. 
Since 3 Na + enter the cell for every Ca 2+ extruded, this 
increase in driving force increases the inward membrane 
current which prolongs the APD. Therefore, depending 
on the relative contributions of lea arl d ^NaCa, increasing 
the magnitude of the calcium transient can either prolong 
(positive coupling) or shorten (negative coupling) the 
APD, as illustrated in Fig. [21 The sign of this coupling 
can be changed in the model by varying the exponent 7 in 
the phenomenological expression f^ a = 1/ [1 + (c s /c s ) 7 ] 
for the steady-state value of /c a , where the constant c s 
sets the concentration range for inactivation. Increasing 
7 enhances calcium-dependent inactivation of 7ca an d 
tends to make the coupling negative. 

Numerical results. The dynamics of the model was 
studied numerically as a function of the two instability 
parameters u and 77 which promote Ca 2+ and voltage 
alternans, respectively, and for two values of 7 that were 
found to yield a positive (7 = 0.7) and a negative (7 = 
1.5) coupling between voltage and calcium. All the other 
parameters are the same as in Ref. 0,^2 an d the pacing 
period is fixed to T = 300 ms. 

The results plotted in Fig. highlight the crucial role 
of the coupling between voltage and calcium in the dy- 
namics. For positive coupling, the instability of the 1:1 
periodic state always occurs through a period-doubling 
bifurcation to clcctromcchanically concordant alternans 
with the long (short) APD corresponding to a large 
(small) calcium transient, independently of whether volt- 
age or calcium is the dominant instability mechanism. In 
contrast, for negative coupling, three distinct modes of 
instability arc found that correspond to (i) concordant 
alternans, as for positive coupling, but only when the 
instability is dominated by voltage (large r/ and small 
u), (ii) discordant alternans with the long (short) APD 
corresponding to a small (large) calcium transient when 
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FIG. 3: Stability boundaries in the ionic model for positive 
(dashed line; 7 = 0.7) and negative (solid line; 7 = 1.5) 
coupling. T = 300 ms. Examples of steady-state dynamics 
close to the stability boundaries are illustrated by plots of 
peak calcium concentration (c^ cak ) vs. APD for a few labelled 
points. Higher order periodicities and irregular dynamics are 
observed further away from these boundaries. 




FIG. 4: Definition of map variables. 
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FIG. 5: Stability boundaries from the map analysis for pos- 
itive coupling C — 0.1 with concordant alternans along the 
dashed line, and negative coupling C = —0.1 (solid line) with 
along the segments a — b: concordant alternans; c — d: discor- 
dant alternans, and 6 — c: quasiperiodicity. 



the instability is dominated by calcium (small Tf and 
large u) , and (iii) quasiperiodic oscillations of APD and 
calcium transient amplitude with a phase and a Hopf 
frequency that vary with Tf and u for the in between 
case where the instability is driven by both voltage and 
calcium. Both electromechanically concordant and dis- 
cordant alternans have been widely observed experimen- 
tally under various conditions |2(ij |. In addition, there is 
experimental evidence for quasiperiodicity in recordings 
of voltage |2l| and, more recently calcium |22^ . 

Iterated map of voltage- calcium dynamics. To inter- 
pret our results, we extend the two-dimensional iterated 
map developed in Ref. for calcium cycling when the 
cell is paced with a fixed periodic voltage waveform, to 
the present case where the voltage is undamped. To a 
good approximation, c s w Cj and cL w Cj preceding a 
stimulus such that we only need to track beat-to- 
beat changes of Cj and Cj. Furthermore, we assume for 
simplicity that buffering of calcium is instantaneous such 
that there exists a unique nonlinear relationship between 
the concentration of free calcium a (cj ) and total calcium 
(free plus bound) cj {cj). The basic variables of the map 
(Fig. f4J are then cf and cj at time t n = nT of the n th + 1 
stimulus, defined by x n = cf(t„) and y n = (v sr /vi)cj (t n ) 
where both x n and y n are in units of /xM, and the APD 
corresponding to this stimulus, A n+ \. 



The map is obtained by extending the restitution map 
to include the effect of calcium on the APD and by inte- 



grating the calcium flux equations 

cj = I T el — Iixp — iGa, + -^NaCa; (9) 

cj = (Vi/V ST )(-I le l + I up ), (10) 

from time t n to time t n+ \. This yields 

A n +i = F(D n ,x n ,y n ), (11) 

x n+ i = x n + R n - U n + A n , (12) 

Vn+l = Vn - Rn + U ni (13) 



respectively, where R n , U n , and A„ are the integrals 
of irei, J up , and — Jca + ^NaCa over the time interval 
[t n ,t n+ i], respectively, and are functions of (D n ,x n ,y n ) 
for a fixed pacing period; ViR n and u,Z7„ are the total 
amount of calcium released from and pumped into the 
SR over one beat, respectively, and ViA n is the net total 
calcium entry into the cell over one beat which can be 
positive (negative) if the exchanger extrudes more (less) 
calcium from the cell than Tea brings into the cell. 

To study the stability of the fixed point of the map 
(A* ,x* ,y*), we exploit the fact that the total amount of 
calcium inside the cell is approximately constant during 
steady-state pacing. Hence, we can approximate the 3- 
dimensional (3-d) map (Eqs. I11I13|) by a 2-d map by 



4 



assuming that z n « z* , where ViZ n = Vi(x n + y n ) is the 
total calcium in the cell at time t n . This 2-d map is 
given by Eqs. HT1 and IHl with D n = T — A n , A„ = 0, and 
y n = z* — x n . A straightforward linear stability analysis 
of this 2-d map yields the eigenvalues 



-A„ -K± ^(K-KY+tC 



where we have defined the quantities 



At, — 
A c = 

C = 



8F/dD n , 

_ x _ d(R n - U n ) d(R n - U n ) 

dx n dy n 

d{R„ - U„) ( OF OF 

9D n \dy n dx n 



(14) 

(15) 
(16) 

(17) 



which are evaluated at the fixed point of the map. Here, 
At, and A c govern the degree of instability of the voltage 
and calcium systems, respectively, while C determines 
the sign of the coupling between the two systems. Mak- 
ing APD-restitution (dF/dD n ) or the relationship be- 
tween release and SR-load (dR n /dy n ) steeper by increas- 
ing Tf and u in the ionic model is equivalent to increas- 
ing A„ and A c , respectively. Graded release implies that 
d(Rn — U n )/dD n is positive for high pacing rates where 
Jca depends on D n , such that the sign of C is governed 
by OF/ dy n — dF/dx n where the latter reflects the effect 
of the magnitude of the calcium transient on APD via 
Jca and JuaCa (Fig. [2J- The periodic fixed point under- 
goes a period doubling bifurcation when |A_| = 1 and a 
Hopf bifurcation for (A u — A c ) 2 + 4C < when the pair of 
complex eigenvalues X± — r e l ^ ±w \ with r = yO^J^~~C 
and tan lu = —AC — (A c — A„) 2 /(A C + At,), crosses the 



unit circle (r = 1). For the latter case, the beat-to-beat 
oscillations of voltage and calcium are modulated with a 
period 2%/u). Examination of the eigenvectors for C < 
reveals that alternans are discordant when A_ is real and 
A c > A„ . We plot in Fig. |5] the corresponding stabil- 
ity boundaries for positive and negative coupling in the 
(A C ,A„) plane which are remarkably isomorphic to the 
stability boundaries obtained by simulations of the ionic 
model in the (u,Tf) plane of Fig. |3] This agreement 
shows that this simple map captures the main robust 
features of the instability of the voltage-calcium system 
observed in the ionic model and experimentally. 

The numerical study of both the ionic model and the 
map in a nonlinear regime reveals the existence of a rich 
dynamical behavior including higher order periodicities 
(3:3, 4:4, etc) as well as transitions to chaos mediated by 
a period-doubling cascade or intermittency depending on 
the parameters. Moreover, this model naturally contains 
memory 0, [2j| due to the slow change of total calcium 
concentration over several beats. Both of these aspects 
will be discussed in more details elsewhere. 

In conclusion, we have outlined the essential three- 
dimensional parameter space that controls dynamic in- 
stability of membrane voltage coupled to calcium cycling 
and we have presented a theoretical framework in which 
to interpret experiments beyond the limitations of the 
one-dimensional restitution relationship. The main axes 
of this parameter space are the degree of instability of the 
voltage and calcium systems, and the sign of the coupling 
between the two systems, which is an important new pa- 
rameter to emerge from this work. These results offer 
new concepts to help identify the mechanisms which un- 
der ly various heart rhythm disorders. This research is 
supported by NIH SCOR P50-HL52319. 
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